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ABSTRACT 



Aims. The BLLac object RGB J0152+017 (z = 0.080) was predicted to be a very high-energy (VHE; > 100 GeV) y-ray source, due to its high 
X-ray and radio fluxes. Our aim is to understand the radiative processes by investigating the observed emission and its production mechanism 
• ^ • using the High Energy Stereoscopic System (H.E.S.S.) experiment. 

J\. Metiiods. We report recent observations of the BLLac source RGB J0152-I-017 made in late October and November 2007 with the H.E.S.S. 
j array consisting of four imaging atmospheric Cherenkov telescopes. Contemporaneous observations were made in X-rays by the Swift and RXTE 
satellites, in the optical band with the ATOM telescope, and in the radio band with the Nangay Radio Telescope. 

Results. A signal of 173 y-ray photons corresponding to a statistical significance of 6.6 cr was found in the data. The energy spectrum of the 
source can be described by a powerlaw with a spectral index of F = 2.95 + 0.36s,ii, ± 0.20sysf The integral flux above 300 GeV corresponds to ~2% 
of the flux of the Crab nebula. The source spectral energy distribution (SED) can be described using a two-component non-thermal synchrotron 
self-Compton (SSC) leptonic model, except in the optical band, which is dominated by a thermal host galaxy component. The parameters that are 
found are very close to those found in similar SSC studies in TeV blazars. 

Conclusions. RGB J0152+017 is discovered as a source of VHE y-rays by H.E.S.S. The location of its synchrotron peak, as derived from the 
SED in Swift data, allows clear classification as a high-frequency-peaked BLLac (HBL). 

Key words, galaxies: BL Lacertae objects: individual: RGB J0152+017 - gamma rays: observations - galaxies: BL Lacertae objects: general - 
galaxies: active 



1. Introduction 

First detected as a radio source dBecker et all Il991h by the 
NRAO Green Bank Telescope and in the Parkes-MIT-NRAO 
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surveys ([Griffith et all 1 1995b RGB J0152-H017 was lat er iden- 
tified as a BL Lac object by iLaurent-Muehleisen et al.l (Il998h . 
who located it at z = 0.080, and was c laimed as an intermediate- 
freque n cy-peaked BL Lac objec t by iLaurent-Muehleisen et al.l 
(Il999h . iBrinkmann et alj (Il997h report the first detection of 
RGBJ0152H-017 in X-rays in the ROSAT-Giesn Bank (RGB) 
sample. The host is an elliptical galaxy with luminosity Mr = 
-24.0 (iNilsson et al.. .2003,) . The source has high radio and 
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X-ray fluxes, making it a good candidate for VHE emission 
dCostamante & GhiseUinill2002l) . motivating its observation by 
the H.E.S.S. experiment. 

The broad-band SED of BLLac objects is typically char- 
acterised by a double-peaked structure, usually attributed to 
synchrotron radiation in the radio-to-X-ray domain and inverse 
Compton scattering in the y-ray dom ain, which is frequentl y 
explained by SSC models (see, e.g., lAharonian et al.l |2005|). 
However, since the flux of BL Lac objects can be highly vari- 
able (e.g. lKrawczvnski et al.ll20d5 ). stationary versions of these 
models are only relevant for contemporaneous multi-wavelength 
observations of a non-flaring state. The contemporaneous ra- 
dio, optical. X-ray, and VHE observations presented here do not 
show any significant variability, and thus enable the first SSC 
modefling of the emission of RGB J0152+017. 



2. H.E.S.S. observations and results 

RGB J0152H-017 was observed by the H.E.S.S. array consisting 
of four imaging atmospheric Che renkov telescopes, loc ated in 
the Khomas Highland, Namibia (lAharonian et alj |2006a). The 
observations were performed from October 30 to November 14, 
2007. The data were taken in wobble mode, where the telescopes 
point in a direction typically at a n offset of 0.5° from the nom- 
inal target position (lAharonian e t al. 2006a). After applying se- 
lection cuts to the data to reject periods affected by poor weather 
conditions and hardware problems, the total live-time used for 
analysis amounts to 14.7 h. The mean zenith angle of the obser- 
vations is 26.9°. 

The data are calibrated according to lAharonian et al.l (|2004|) . 
Energies are reconstruct ed taking the effective o ptical efficiency 
evolution into account (lAharonian et al.l l2006 a). The separa- 
tion of y-ray-like events from cosmic-ray background events 
was m ade using the Hillas moment-analysis technique dHillasI 
Il985h. Signal extractio n was performed using standard cuts 
(lAharonian et aljr2006al) . The on-source events were taken from 
a circular region around the target with a radius of = 
0.11°. The background was estimated using reflected regions 
(lAharonian et alj |2006a) located at the same offset from the cen- 
tre of the observed field as the on-source region. 

A signal of 173 y-ray events is found from the direction 
of RGB J0152-I-017. The statistical significance of the detection 
is 6.6 cr according to Li & Ma (1983). The preliminary detec- 
tion was reported by iNedbal et alj ( I2007l) . A two-dimensional 
Gaussian fit of the excess yields a position a 72000 = 1^52™33^!5 + 
5^.3sta, ± P.3,sy,t, (Jy20oo = 1°46'40'.'3 + 107^;,, ± 20;'3,,,. The 
measured position is compatible with the nominal position of 
RGBJ0152-H017 (ffyaooo = l''52'"39^78, (Jy20oo = 1°47'18'.'70) 
at the Icr level. Given this spatial coincidence, we identify the 
source of y-rays with RGB J0152-I-017. The angular distribu- 
tion of events coming from RGB J0152H-017, shown in Fig.[Il 
is compatible with the expectation from the Monte Carlo simu- 
lations of a point source. 

Figure|2|shows the time-averaged differential spectrum. The 
spectrum was derived using standard cuts with an energy thresh- 
old of 30 GeV. Another set of cuts, the spectrum cuts de- 
scribed in lAharonian et alj (l2006bl) . were used to lower the en- 
ergy threshold and improve the photon statistics (factor ~2 in- 
crease above the standard cuts). Both give consistent results (see 
inlay in Fig.|2]and caption). Because of the better statistics and 
energy range, we use the spectrum derived using spectrum cuts 
in the following. Between the threshold of 240 GeV and 3.8 TeV, 
this differential spectrum is described well (Y^/d.o.f.=2.16/4) 
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Fig. 1. Angular distribution of excess events. The dot-dashed line 
shows the angular distance cut used for extracting the signal. The 
excess distribution is consistent with the H.E.S.S. point spread 
function as derived from Monte Carlo simulations (solid line). 
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Fig. 2. Differential spectrum of RGB J0152+017. The spectrum 
obtained using spectrum cuts (black closed circles) is compared 
with the one obtained by the standard cuts (blue open circles). 
The black line shows the best fit by a powerlaw function of 
the former. The three points with the highest photon energy 
represent upper limits at 99% confidence level, calculated us- 
ing Feldman & Cousins (1998). All error bars are only statis- 
tical. The fit parameters of a powerlaw fit are F = 2.95 + 
0.36stat ± 0.20syst and (D(lTeV) = (5.7 + 1.6stat ± l-lsyst) x lO^'^ 
cm"^ s ' TeV ' for the spectrum cuts, and F = 3.53 + 0.60stat ± 
0.2syst and <l)(lTeV) = (4.4 + 2.0) x lO"'^ cm^^ s"' TeV"' for the 
standard cuts. The insert shows 1 and 2 cr confidence levels of 
the fit parameters. 



by a power law dN/dE - ^oiE/lTeV) ^ with a photon in- 
dex F - 2.95 + 0.36stat ± 0.20syst and normalisation at 1 TeV 
of a)(lTeV) = (5.7 + 1.6stat + l-Uyst) X lO'^^ cm^^s-i TeV"'. 
The 99% confidence level upper limits f or the highest three 
bins sh own in Fig. |2]were calculated using iFeldman & CousinsI 
(119981) . 

The integral flux above 300 GeV is / = (2.70 + 0.51stat ± 
0.54syst) X 10^'^cm"^ s"', which corresponds to -2% of the flux 
of the Crab nebula abov e the same threshold as determined by 
lAharonian et al.l(l2006ah . Figure|3]shows the nightly evolution of 
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Fig. 3. Mean nightly integral flux from RGB JO 152 -1-017 above 
300 GeV. Only the statistical errors are shown. Upper limits at 
99% confidence level are calculated when no signal is found 
(grey points). The dashed line shows a fit of a constant to the 
data points with;^f^/d.o.f. of 17.2/12. The fit was performed us- 
ing all nights. 

the y-ray flux above 300 GeV. There is no significant variability 
between nights in the lightcurve. The ;^^^/d.o.f. of the fit to a 
constant is 17.2/12, corresponding to ax^ probability of 14%. 

All results were checked with independent analysis proce- 
dures and calibration chain giving consistent results. 

3. Multi-wavelength observations with Swift, RXTE, 
ATOM, and the Nangay Radio Telescope 

3.1. X-ray data from Swift and RXTE 

Target of opportunity (ToO) observations of RGB J0152+017 
were performed with Swift and RXTE on November 13, 14, and 
15, 2007 ti-iggered by the H.E.S.S. discovery. 

The Swift/XRT (Buitows et al. 2005) data (5.44ks) were 
taken in photon-counting mode. The spectra were extracted with 
xselect v2 . 4 from a circular region with a radius of 20 pixels 
(0.8') around the position of RGB JO 152-1-0 17, which contains 
90% of the PSF at 1 .5 keV. An appropriate background was ex- 
tracted from a region next to the source with four times this 
area. The auxiliary response files were created with the script 
xrtmkarf vQ . 5 . 6 and the response matrices were taken from 
the Swift package of the calibration database caldb v3.4.1. 
Due to the low count rate of 0.3 cts/s, any pileup effect on the 
spectrum is negligible. We find no significant variability dur- 
ing any of the pointings or between the three subsequent days; 
hence, individual spectra were combined to achieve better pho- 
ton statistics. The spectral analysis was performed using the tool 
Xspec vll.3.2. Abrokenpowerlaw (Ti = 1.93 + 0.20, r2 = 
2.82 + 0.13,£'bi-eak = 1-29 + 0.12keV) wi th a Galactic absorp- 
tion of 2.72 X lO^^cm ^ (LAB Survey. iKalberla et alj llOOSh 
is a good description (x^/d.o.f. = 24/26), and the resulting 
unabsorbed flux is /^().5-2keV ~ 5.1 x 10"^^ erg cm"^ s"' and 
■F2-iokeV ~ 2.7 X 10^'^ergcm^^s"^ 

Simultaneous observations at hig her X-ray photon e nergies 
were obtained with the RXTE/PCA dJahoda et al.lll996h . Only 
data of PCU2 and the top layer were taken to obtain the best 
signal-to-noise ratio. After filtering out the influence of the South 
Atlantic Anomaly, tracking offsets, and the electron contamina- 
tion, an exposure of 3.2 ks remains. Given the low count rate of 
0.7 cts/s, the "faint background model" provided by the RXTE 
Guest Observer Facility was used to generate the background 
spectrum with the script pcabackest v3 . 1. The response ma- 
trices were created with pcarsp vlO.l. Again no significant 
variations were found between the three observations, and indi- 



vidual spectra were combined to achieve better photon statistics. 
The PCA spectrum can be described by an absorbed single pow- 
erlaw with photon index F = 2.72 + 0.08 (x'-/d.o.f. = 20/16) 
between 2 and lOkeV, using the same Galactic absorption as for 
Swift data. The resulting flux F2-iokeV ~ 6.8 x 10"'^ erg cm"^ s"' 
exceeds the one obtained simultaneously with Swift by a factor 
of 2.5. We attribute this mostly to contamination by the nearby 
galaxy cluster Abell267 (44.6' offset from RGB J0152+017 but 
still in the field of view of the non-imaging PCA). 

A detailed decomposition is beyond the scope of this paper, 
so we exclude RXTE data from broadband modelling. The RXTE 
data-set confirms the absence of variability during November 
2007, also in the energy range up to 10 keV. For the SED mod- 
elling, the average spectrum is treated as an upper limit. Further 
observations with RXTE in December 2007 also show no indica- 
tion of variability. 

3.2. Optical data 

Optical observations were taken using the ATOM telescope 
( Hauser et al. 2004 ) at the H.E.S.S site from November 10, 2007. 
No significant variability was detected during the nights between 
November 10 and November 20; R-band fluxes binned nightly 
show an RMS of 0.02 mag. 

Absolute flux values were found using differential photome- 
try against stars calibrated by K. Nilsson (priv. comm.). We mea- 
sured a total flux of = 15.25 + 0.01 mag (host galaxy + core) 
in an aperture of 4" radius. The host galaxy was subtracted with 
galaxy parameters given in Nilsson et al. (2003), and aperture 
correction given in Eq. (4) of lYound (119761) . The core flux in the 
R-band (640 nm) was found to be 0.62 + 0.08 mJy. This value 
was not corrected for Galactic extinction. 

3.3. Radio data 

The Nangay Radio Telescope (NRT) is a meridian transit 
telescope with a m ain spherical mirror of 300 m x 35 m 
dTheureau et al.l20^7l) . The low-frequency receiver, covering the 
band 1.8-3.5 GHz was used, with the NRT standard filterbank 
backend. 

The NRT observations were obtained in two contiguous 
bands of 12.5 MHz bandwidth, centred at 2679 and 2691 MHz 
(average frequency: 2685 MHz). Two linear polarisation re- 
ceivers were used during the 22 60-second drift scan observa- 
tions on the source on November 12 and 14, 2007. The data have 
been processed with the standard NRT software packages NAPS 
and SIR. All bands and polarisations have been averaged, giving 
an RMS noise of 2.2 mJy. The source 3C 295 was observed for 
calibration, on November 11, 13, and 15, 2007. 

Taking into account a flux density for this source of 1 2.30 ± 
0.06 Jy using the spectral fit published by lOtt et all d 19941) . 
we derived a flux density of 56 + 6mJy at 2685 MHz for 
RGB J0152H-017. No significant variability was found in the ra- 
dio data. 



4. Discussion 

Figure H shows the SED of RGB J0152-H017 with the data from 
Nangay, ATOM, Swift/XRT, RXTE/PCA, and H.E.S.S. Even 
though some data are not strictly simultaneous, no significant 
variability is found in the X-ray and optical bands throughout 
the periods covered; hence, a common modelling of the contem- 
poraneous X-ray and VHE data appears justified. 
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Fig. 4. The spectral energy distribution of 
RGBJ0152+017. Shown are the H.E.S.S. 
spectrum (red filled circles and upper lim- 
its), and contemporaneous RXTE (blue open 
triangles), SwiftlXKI (corrected for Galactic 
absorption, magenta filled circles), opti- 
cal host galaxy-subtracted (ATOM) and ra- 
dio (Nangay) observations (large red filled 
squares). The black crosses are archival data. 
The blue open points in the optical R-band 
corre spond to the total a nd the core fluxes 
from lNilsson etalj (l2003h . A blob -m-jet syn- 
chrotron self-Compton model (see text) ap- 
pHed to RGBJ0152+017 is also shown, de- 
scribing the soft X-ray and VHE parts of the 
SED, with a simple synchrotron model shown 
at low frequencies to describe the extended 
part of the jet. The contribution of the dom- 
inating host galaxy is shown in the optical 
band. The dashed line above the solid line at 
VHE shows the source spectrum after correct- 
ing for EBL absorption. The left- and right- 
hand side inlays detail portions of the ob- 
served X-ray and VHE specti-um, respectively. 



The optical part of the SED is mainly due to the host 
galaxy, which is de tected and resolved in optical wavelengths 
jNilsson et al.ll2003l) . A template of the spectrum of such an el- 
liptic galaxy is shown in the SED, as infe rred from the code 
PEGASE (.Fioc & Rocca-Volmerang3ll997h . The host-galaxy- 
subtracted data point from the ATOM telescope might include 
several additional contributions — from an accretion disk, an ex- 
tended jet (see below), or a central stellar population — so that it 
is considered as an upper limit in the following SSC model. A 
model including the optical ATOM data with possible additional 
contributions is beyond the scope of this paper. 

We applied a non-thermal leptonic SSC model 

dKatarzvriski et al.l 1200 Ih to account for the contemporane- 
ous observations by Swift in X-rays and by H.E.S.S. in the 
VHE band. The radio data are assumed to originate in an 
extended region, described by a separate synchrotron model 
for the extended jet (iKatarzyiiski et alj l200l1) to explain the 
low-frequen cy part of the SED (as in, e.g., lAharonian et al.l 
I2005ll2008l) . 

We should emphasise that the aim of applying this model 
in this work is not to present a definitive interpretation for this 
source, but rather to show that a standard SSC model is able to 
account for the VHE and Swift X-ray observations. 

For the SSC model, we describe the system as a small ho- 
mogeneous spherical, emitting region (blob) of radius R within 
the jet, filled with a tangled magnetic field B and propagating 
with a Doppler factor 5 - [r(l -yScosfl)] Here F is the bulk 
Lorentz factor of the emitting plasma blob, (3 = v/c, and 9 is 
the angle of the velocity vector, with respect to the line-of-sight. 
The electron energy distribution (FED) is described by a broken 
powerlaw, with indices n\ and n2, between Lorentz factors y^in 
and jmux, with a break at ybieak and density normalisation K. 

The model also accounts for the absorption by the extra- 
galactic background light (EBL) with the parameters given in 
IPrimack et al. (2005). RGB J0152-H017 is too nearby (z = 0.08) 
to add to the constraints on the E BL that were found by H .E.S.S. 
measurements of other blazars ( Aharonian et alj|2006d) . In all 



the models, we assume Hq = 70kms-iMpc-i, giving a lumi- 
nosity distance of t/z, = 1.078 x 10^^ cm for RGB J0152+017. 

The FED can be described by = 3.1 x 10'*cm 
Tmin = 1, and ymax =4 X lO*". The break energy is assumed 
at ybreak = 7.0 X 10^ and is consistent with the Swift [XKl spec- 
trum, while providing good agreement with the H.E.S.S. data. 
We assume the canonical index n\ =2.0 for the low-energy part 
of the FED, in accordance with standard Fermi-type acceleration 
mechanisms. The value «2 = 3.0 for the high-energy part of the 
FED is constrained by the high-energy part of the X-ray spec- 
trum. A good solution is found with the emitting region charac- 
terised by 5 = 25, 1.5 X lO'^cm, andB = O.IOG. 

For the extended jet, the data are described well by R^^i - 
lO'^cm, 5jet = 7, /iTjet = 70cm"^ Bjet = 0.05 G, and ybieakjet = 
10"* at the bas e of the jet, and Lfet — 50 pc (all the parameters are 
detailed in Kat arzyriski et al.| |200lh. 

Assuming additional contributions in the optical band, the 
multi-wavelength SED can thus be explained with a standard 
shock-acceleration process. The parameters derived from the 
mod el are similar to previo us results for this type of source (see, 
e.g.. lGhisellini et al.ll2002l) . 

From the current Nan9ay radio data and the Swift X-ray data, 
we obtain a broad-band spectral index a„ ~ 0.56 between the 
radio and the X-ray domains. The obtained SED, the correspond- 
ing location of the synchrotron peak, and the flux and shape of 
the Swift spectrum lead us to conclude that RGB J0152H-017 can 
clearly be classified as an HBL object at the time of H.E.S.S. 
observations. 



5. Conclusion 

The HBL RGB J0152H-017 was detected in VHE at energies 
> 300 GeV with the H.E.S.S. experiment. The contemporane- 
ous Swift, RXTE, Nan^ay, ATOM, and H.E.S.S. data allow the 
multi-wavelength SED for RGB J0152-I-017 to be derived for 
the first time , and to clearly confirm its HBL nature at the time 
of the H.E.S.S. observations. In general, large variations of the 
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VHE flux are expected in TeV blazars, making further monitor- 
ing of this source to detect high states of the VHE flux (flares) 
desirable. 
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